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Abstract. We consider some nonlinear phenomena in metamaterials with nega- 
tive refractive index properties. Our consideration includes a survey of previously 
known results as well as identification of the phenomena that are important for 
applications of this new field. We focus on optical behavior of thin films as well 
as multi-wave interactions. 

1 Introduction 

In recent years, the development of nanotechnology has led to the creation of new materials 
with very unusual optical properties - metamaterials. Metamaterials are composed of host di- 
electric materials with embedded periodic structures, patterned on a scale much shorter than 
the operating wavelength and thus on a nanometer scale if the infrared or visible spectra are 
targeted. Though they are made of positive index materials at small length scale, these struc- 
tures exhibit abnormal dispersion characteristics at optical scales. The material can behave 
as a medium with negative index of refraction associated with negative values of the effective 
permittivity e and effective permeability fj, p]. The electric field, magnetic field wavevector then 
constitute a left handed system, where the wavevector and the Poynting vector are antiparal- 
lel pf3 4 5 . In this case, they are termed left handed materials (LHMs). Beyond this, inversion 
of the phase velocity with respect to the direction of propagation of energy leads to negative 
refraction. That means both the incident and refracted beam are on the same side of the normal 
to the refraction interface, as if predicted by Snell's law with a negative refractive index. 

Left handed properties of metamaterials are most pronounced in two situations. The first 
is in the behavior of a quasi-monochromatic wave at the interface of a left-handed and a 
conventional right handed material. In particular we consider layered materials because the 
use of thin films of LHM opens new opportunities for design of photonic crystals with unusual 
properties. The second is the interaction of two waves with different carrier frequencies which 
experience different signs of refractive index. Such a difference is a consequence of how LHM are 
realized using simultaneous resonance of electric and magnetic field components with metallic 
nanostructures embedded in a host dielectric. Such a material behaves as right handed for waves 
with carrier frequencies above the resonance region and as left-handed for waves with carrier 
frequencies below the resonance region. 

In this paper we consider only phenomena related to these two cases. 
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2 Nonlinear waveguiding phenomena on the interface of materials with 
opposite signs of their refractive indices 

Let us consider a layer of left-handed material on a dielectric substrate with positive refraction 
index - a right handed material (RHM) . We study optical waves propagating along this interface 
in the z-direction of the Cartesian frame. The normal to the surface is chosen to be the x- 
direction while the y-axis lies in the interface plane. There are two types of the surface waves: 
(a) transverse electric waves (TE waves) with E = {Q.Ey^Q) and H = {Hx,0, Hz), and (b) 
transverse magnetic waves (TM waves) with E = {Ex,0,Ez) and H — {O.Hy,0). If the tensor 
of dielectric permeability of one of the media is a square-law function of electrical field strength 
e{uj, E) = e{uj) + e„i(a;) : EE then propagation of the non-hnear surface polaritons (NLSP) is 
possible for both the TM and the TE waves |10I12I14| . Linear surface polaritons corresponding 
to e„i(a^) — were found in f7^. Without loss of generality, here we will consider only TE waves. 



2.1 Non-linear surface waves of TE-type 

Let us consider waves on the interface between two isotropic media. The electrodynamic prop- 
erties of each medium are described by effective dielectric permittivity and effective magnetic 
permeability. In the linear limit at x < the medium is characterized by positive permittivity 
and permeability ei > 0,/Lti > and at a; > both e and /i2 < are chosen to be negative 
62 < 0, /i2 < 0. Due to its nature the electrical field vector of the TE-waves has only one 
nonzero component, Ey — E, which is governed by the equation: 

E^zz+E^^, + kle{LJ,-E)fi{uj,-E)E = 0, (1) 

where fcg = uj/c. The magnetic field vector components are coupled with E by the following 
relations: 

H., = i{ko^i)-^E^z, H, = -i{kn^i)-^E^, 

On the interface plane at a; = the field components E, and satisfy the set of jump 
conditions: 

E{0-) = E{0+), HziO~) = HziO+), /iii/,(0-) = M2i/.(0+). 

Due to the translation symmetry along the z-axis, the y-component of the electric field E{u;, x, z) 
can be chosen as 

E{u), X, z) = <P{uj, x) exp[i/3(w)2:], (2) 

where /3(w) is a propagation constant [5]. The equation for the transverse profile of the electrical 
field of NLSW follows from equations dJ) and ©: 

^^xa= + [kle{Lu,x)ii{uj,x)- p'']^ = 0. (3) 

The above jump conditions now can be represented in terms of the amplitude function ^(x): 

<^(0-) = <Z.(0+), /ir'^,x(0-)=M2"^<?,x(0+), (4) 

The boundary conditions 

lim = 0, lim ^^^(tj,x) = (5) 

— >oo \x\ — >oo 

select the surface wave solutions of Q. 

As an example giving rise to nonlinear polaritons we considered the simple geometry where 
the nonlinear left-handed material is placed on a linear semi-infinite conventional right-handed 
material: 

/ N / £i, X <Q , . / ^1, x < 

^ 1 £2 + Snl'P , X > ' ' ^ \fi2,X>0 
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Equation ([3]) can be solved at a; < and at a; > 0, where the media are homogeneous, taking 
account of the boundary conditions ([5]). The matching condition for these solutions at the 
boundary a; = results in dispersion relations for NLSP. 

The equation ([3]) does not have bounded solutions, if fi2eni < (self-refocusing LHM). This 
means that nonlinear surface polaritons do not exist in this geometry with cubic nonlinearity in 
presence of a self-refocusing left-handed material. It can be shown that, if /i2£ni > (this is the 
case of a self- focusing non-linear LHM), then NLSP exists only when = (/3^ — fcoe2M2) > 0. 
In this case the solution of equation ^ reads 

^ r Aiexp(gx), a; < 

[a^'^/^ scch[p{x - X2)], X > 0, ^ ' 

where a — fcge2M2/2p^, 9^ = [0^ — /coEi/ii) > 0, and X2 is the location of the maximum of the 
electrical field. The continuity condition at the interface x ~ Q requires 

^1 = a"^/^ sech(px2), (7) 

l^2q = MiP tanh(pa;2). (8) 

The expression (O defines X2 as function of electric field strength Ai at the interface. This 
equation ([8]) is the dispersion relation for NLSP. Since mui and ^2 have different signs ( /ii > 0, 
M2 < 0), the parameter X2 must be negative. Thus the electric field reaches its maximum value 
on the interface at a: = 0. It should be noted that in the case where both linear and nonlinear 
dielectrics are right-handed, the electric field reaches its maximum value inside the non-linear 
dielectric medium, i.e., where X2 is positive. 

In the other case, when a linear LHM is placed on a nonlinear semi-infinite medium (RHM) , 
e and /i are defined as follows: 

82, x>0' ^^^'^)=|/i2,x>0 

The governing equations and corresponding solutions in this case can be easily obtained from 
equations ® and d?]), dH]). 

For a self-refocusing medium, solutions of equation ([3]) are given by the expression 

<^{x) = [ ±l"r^^^cosech[g(a; - xi)], x<0 ^^-^ 
I exp( — p^), X ^ 

The continuity condition at the interface x = Q requires 

A2 ^ ±\a\-^''^ I aiiM-qxi) (10) 

^ig = ^2'7tanh(ga;i) (11) 

Since /.ti > and /i2 < 0, it follows from the dispersion relation pip that the parameter xi 
must be negative. Thus the electric field reaches its maximum value on the interface a; = as 
in the case of a self-refocusing medium. 

A detailed analysis of the surface waves on the interface between nonlinear LHM and non- 
linear RHM was presented in [14|llj . 



2.2 Nonlinear planar waveguide 

In this section we consider a sandwich type structure, which consists of a dielectric layer with 
left-handed properties surrounded by Kerr-like nonlinear dielectrics |13) . The presence of two 
interfaces in this configuration (at a; = and x — h) could lead to two types of guided waves. 
The first type is represented by an internal wave localized in the layer. We will refer to these 
types of solutions as non- linear guided modes (NLGM). The second type is made up of two 
coupled surface waves and they will be referred to as non- linear surface waves (NLSW) . 
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2.2.1 Dispersion relations of nonlinear surface TE-waves 

The equation ([3]) describing the transverse electric field profile of the wave in a three-layer 
structure can be decomposed into a system of simple equations. If the condition ^.isfj > 
holds, then the solution of these equations is 

{Qf]^ sech[pi(a; — xi)], x <Q 

^exp(— Kx) + i3exp(Kx), Q <x <h , 
a2 sech[p2(x — X2)\, x>h 

where — f3'^ — /cqEi^i, p^ = (3'^ — kQ£2H2, i — 1,2 and n'^ — — klssfJ-s > 0. The signs oi pf 
and P2 are chosen to be positive to guarantee validity of the boundary conditions ([5]). The sign 
of is also chosen to be positive (see (O), as we will see below such a choice corresponds to 

the case of surface waves. For ^J-ie^^l < we have 

{±|Q;i|^-^/^cosech[pi(a:: — xi)], x<0 
Aexp{—Kx) + Bexp{Kx), 0<x<h . 
±|a2|~"^^^cosech[p2(a:^ — a;2)], x>h 

There xi^2 are the constants of integration corresponding to the co-ordinates of the maximum of 
the transverse profile of the electric field. The continuity conditions of the electric and magnetic 
field components at the interfaces ([4]) result in a homogeneous system of linear equations. There 
are nontrivial solutions of this system if the determinant is zero, which leads to the following 
dispersion relation: 

e^"'* (1 + AisPi/K) (1 + MsP2/k) = (1 - (J-sPi/k) (1 - ^XsP2/k) , (12) 

where pi,2 — Ai]^2Pi.2 tanh(pi^22;i,2)- Introducing new variables (f)i and 02 as tanh(0i/2) = 
fisPi/K, tanh(02/2) = Ms?'2/'*i we represent relation p2p in the form: 

2k/i + 01 + 02 = 0. (13) 

Since nh is positive, at least one of 0i and 02 must be negative. We define the amplitudes 

— 1/2 

of nonlinear surface waves at the corresponding interfaces as: Agi = sech(pia;i), As2 — 

1/2 

a2 sech[p2(^ — 2:2)]- The connection between parameters Xi and X2 can be found from the 
condition describing coupling of these two surface waves: 

As2 = Asi [cosh Kh — {^JLsP2l 1^) smYinh] . (14) 

Taking into account the relation (fT4| we conclude that the dispersion relation ([12]) contains xi 
as a free parameter. On the other hand, parameter Asi is defined in terms of xi, therefore the 
value of the electric field at a; = can be chosen as a free parameter instead of xi. Finally, 
equation (|12|) can be considered as an implicit one parameter expression for describing the 
dependance of the propagation constant on the frequency f3 = /3(a;; Agi). 



2.2.2 Dispersion relations of nonlinear guided waves 

Analysis of the nonlinear guided wave is based on equation In contrast to the previous case, 
here we choose the sign of the coefficient in front of the function <P in the equation ([3]) to be 
negative i.e. ^qEs/Zs — = > 0. If fJ-isllj > 0, the solution of equation ^ is represented by 
the expression 

{ai^^^ sechlpi^x — xi)], x<0 
Aeyip{iKix) + S exp(— zkix), < x < h . 
— 1/2 
a2 sech\p2{x — X2)], x>h 
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or 



In this case the dispersion relation for NLGW reads as 

(1 + ifisPi/ni) (1 + 'if^sP2/Ki) 

tan(Kih) = '1 . (16) 

If we define the phases 4>i and (j)2 by the formulae: tan((/)i/2) = /i^pi/Ki, tan((/)2/2) = /isi52/'*i, 
then equation ([2]) can be written as 

2Kih = 4>i+(f>2 + 2TTm, TO = 0, 1, 2, ... (17) 

As in to the case of linear waveguides, this expression shows that the full phase shift of the 
zigzag wave |8| consists of contributions from the linear medium of the LHM dielectric slab 
and the phase shifts — 0i and —02, which occur in the total internal reflection at the linear- 
non-linear interfaces. This nonlinear phase shift can be named the non-linear Goos-Hanchen 
effect |17I18I19| 

Analysis of TM waves is presented in [10115] . The corresponding dispersion relation and 
NLGW modes were studied in [14j . In particular, symmetric, asymmetric, and antisymmetric, 
forward and backward modes have been found. 



2.3 Wave refraction and reflection at tliin-film on interface between two dielectrics 

Let us consider the case where a thin film is inserted between two isotropic media. The optical 
properties of these media are described by the permittivity and permeability: £i and /xi at 
X < , 62 and /i2 at x > 0. Let us assume that the thin film of metamaterial separates 
these two media. The optical properties of the metamaterial are described by the polarization 
P'^''-' and the magnetization M^*'. We assume that the width of film If is much less than the 
carrier wave length of a quasiharmonic electromagnetic wave. To describe the propagation of the 
electromagnetic wave through thin film it is sufficient to take account of the jump in the electric 
and magnetic field components induced by polarization and magnetization of the film. Following 
to |20I21| we can approximate transverse distribution of polarization and magnetization in the 
film by a delta-function located at a; = 0. This approximation leads to the following jump 
conditions 

EyiO~)~EyiO+) = ^^^M(/\ i7,(0-)-i/,(0+) = ^|pi^), (18) 
for the TE wave, and 

F,(0-) - Hy{0+) = -^§-^P^'\ E^iO-) - E^{0+) = -^^Mi^\ (19a) 

for the TM wave. The expressions ^5)1 and p9ap are modifications the corresponding bound- 
ary relations used in ' 20121122123124] . These new relations are applicable both for the case of 
continuum waves and for ultra-short or extremely short pulses. 

It follows from and (|19ap that the electric and magnetic field components have different 
values in the media surrounding the film. The values of E and H (i.e., field values acting on 
nanostructures (meta-atoms) of the film), are naturally defined as follows: 

Ea{x = 0) = 0.5 [Ea{{)-) + Ea{0~)] (20) 

Ha{x = 0) = 0.5 [-ffa(O-) + i?a(0-)] (21) 

where a = y or z. In the case of a nonmagnetic film, using (|18p . we obtain the well known 
result pin^ for a TE-wave Ey{x = 0) = ^^^(0-) = ^^^(0-). 
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Let us consider a plane wave which has normal incidence to the interface. For the sake 
of simplicity we assume that the two dielectric media surrounding the film are dispersionless. 
Then in the general case the electric field can be represented in terms of incident, reflected and 
transmitted waves: 

E(x t)^( ^"(^ " ^/^i^ + ^'■^/(^ + ^/^i^ ^ < ^ (22) 
^ ' ' \ Etr{t- X/V2) a; > ' ^ ' 

Here Vi is group velocity in the dielectric at a; < 0, and V2 is group velocity in dielectric at 
X > 0. Using the jump conditions (|19ap and equation Ey^^ = —c~^Hz^t we find the analog 

of the Fresnel relations: 

E At) - ^^E- + i^M(^) - :^^P(4 (24) 

where Py*'' = Py{t, x = 0) and Afi^"* = Mz{t, x — 0) are surface polarization and magnetization 
of the thin film respectively. 

For the TM-wave, using a similar approach, we can obtain following the Fresnel relations: 

eiVi+e2V2 c(£i\/i + £2^2) Lot ^ e iVi ot J 

where the magnetic field amplitudes are defined by: 

„ ( ,N / H^t - x/Vi) + Href{t + x/Fi), a; < 

^^^"^'^^"1 i?t,(t-x/l/2), X>0- 

The electric field can be defined from the equation Hy ,^ = c^^E^.t- 

In what follows we will consider the case of a TE wave since there is no differences between 
TE and TM waves in case of normal incident waves. To obtain the polarization and magne- 
tization of the thin film it is necessary to choose an appropriate model for the meta-atoms of 
the film. We start from the hypotheses that the linear response of film is determined by the 
following expressions: e/(w) = 1 — {ajp/Lu)"^, fJ.f{uj) = 1 — /Sm^^/ {uj^ — uj'^)- This response follows 
from a standard linear Lorenz model which considers meta-atoms as linear oscillators forced by 
electric and magnetic fields. One natural generalization of such a model is to take into account 
anharmonicity of the electrical oscillations in nanostructures. This leads to the simplest model 
for nonlinear response of metamaterials: 

Ptt + uj^P + TePt + gpP3 = {ujl/An)E, (27) 

M^tt + UjIP + PrrM^ = -(/3™/4^)i/,tt, (28) 

where cud is frequency of the dimensional quantization due to confinement of the plasma in the 
nanostructures. Losses in the metallic nanostructures are taken into account by the parameters 
Pe and Pm . Let us normalize the variables in equations ((23l [24)1 and ((27l [28]) 

etr = Etr/^TrPo,eref = £^re/ /47rPo , ei„ = Pi„/47rPo, 

q = P/Po, m = M/(3mV^Po, T = ujpt. 
The system of the normalized equations takes the following form 

/ ^ IfiOv ( 9q ^ , dm'] 
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efir) 



IfLU, 



fUJp 



dq 



dm 



m,rT + iyiTl^p)^ra + {r,n/ U0p)m^r 



-0.5 



etr(T) + \Je\le2 (em(T) + e,.e/(T)) 



(30) 



(31) 
(32) 



Inhomogcncous broadening of the resonance line can be taken into account if we replace q and m 
in by (g) and (to). Here corner brackets (••) denote averaging over different frequencies 

iOd and LOT- 

We assume that the pulse duration is much shorter than the characteristic plasmonic oscil- 
lation damping time, and neglect the dissipation terms in (j29l - I32p . Furthermore, if the slowly 
varying envelope approximation is acceptable, we can reduce the system of these equations and 
consider the following system 



£tr{t) 



2qi{ujQ) 



qi{uJo) + 92(^0) 



f. 



o /,N _ g2(^o)-gl(c^o) ^ 



qi{uJo) + q2{uJo) 



[ko {V)+qiiLUo){M)], 

[ko {V) - q2{uJo) {M)] 



gi(tJo) + 92(^^0) '?i('^o) + 92(^0) 



(33) 
(34) 

(35) 
(36) 



Here fco = ^o/c, loq is the carrier wave frequency, Sin, £tr and £ref are envelopes of the incident, 
transmitted and reflected pulses, V and JV[ are envelopes of the polarization and magnetization. 
The wavenumbers gi,2 are determined as gi,2 = ko^ei^ cos ^1,2; Qi and O2 are incident and 
refraction angles respectively; Uat is volume density of meta-atoms. The envelopes of electric 
£f and magnetic 7i/ fields, acting on nanostructures (meta-atoms) of the thin film can be 
represented as follows: 

£f — [£tr + £in + ^ref] /2, H f — \\fe2£tr + \f^{£in ~ £ref)\ /2. 

Using the normalized variables etr = ^tr/^O: e^e/ — £ref /Aq, e-m — £in/Ao, q = V / Po^m = 
M./Mq,t = iOpt, and introducing parameters Pq = (wp/87rwo)^o andAfo — (/3mWoY^/87rcjp)Ao, 
we present the dimensionless system of equations describing the interaction of the electromag- 
netic pulse with a thin film of metamaterial 



etr{t) 
Crefit) 



2ni 



ikoUatl 



f 



ni + 712 
n2 - ni 
ni + n2 



2(ni -I- 712) 
2(ni -I- 712) 



— (9) + /3m"l7l2— (m) 



P 

LOO 

LO- 



LOv 



- (9) - /?m«2»^2— {m) 



OJQ 



(37) 
(38) 



q.T = ii^p ^(t^o - LUd)q + ig2\q?q + i {etr + ein + ere/) /2, (39) 
m^r = iujp^iujQ - cjt)to + i [etr + (ni/"2)(em - eref)] /2. (40) 

Here 711.2 = ^£1,2 are refractive indices of the dielectric media. The constant of anharmonicity 

g2, used in these equations, is defined as 52 = (3(7p/2a;oWp)(a;p/87ra;o)^^o- The free parameter 
Aq can be represented in terms of the peak intensity of the incident pulse. 

In the case of a thin film with no magnetization the jump conditions (jl8|) results in the 
relation etr = Sin + eref- Equations (|37p -(|40 p in this case transform into the system of equations 
describing the refraction of the ultra-short electromagnetic pulse on the thin film with embedded 
nanoparticles, which was considered in [25j . 
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2.4 Optical bistability in nonlinear layered structures 

A linear layer of dielectric material acts as a resonator for an incident light beam. The transmis- 
sivity of this layer is determined by the thickness and refractive index of the layer material and 
by the wavelength of the incident light. If such a dielectric is nonlinear then its optical thickness 
depends on the intensity of the light field. Therefore, transmissivity of this layer depends on the 
light intensity. Let us consider the resonance reflection condition in the linear limit. Gradual 
increase of light intensity will be followed by a drastic increase in the transmissivity if the light 
intensity reaches a certain threshold level. In other words, the system switches from an opaque 
to a transparent state. This is the phenomena of optical bistability [26]. Phenomena of bistabil- 
ity in thin films were studied in 22 27j . The more realistic case of a film with nonzero thickness 
was considered in |29|28j . The switching condition is determined by the phase difference of the 
waves reflected from the front and back boundaries of the layer. Let us consider a sandwich 
structure with a metamaterial film and a layer of nonlinear dielectric. The presence of the 
metamaterial film results in a change in the effective thickness of the resonator. If the refractive 
index of such a film is negative, then the effective optical thickness of the sandwich-resonator is 
less than the optical thickness of the resonator without film. Therefore, the switching condition 
is realized at a higher level of light intensity [3IJ' . 

It was also shown in [30j that switching light intensity depends on the incident beam direc- 
tion. Switching intensity is higher for the case where the light is entering the structure from 
the side covered by the metamaterial film than for the case of opposite beam direction. 

The hysteresis loop is an important characteristic of bistable devices. Switching thresholds 
from opaque to transparent state and from transparent to opaque state are different. The hys- 
teresis width is the difference of these two switching thresholds. The presence of a metamaterial 
film, as was shown in [3D] , results in an increase of the hysteresis width. 

3 Parametric interactions 

Transformation of frequency of an electromagnetic wave propagating in a non-linear medium 
is one of the most fundamental effects forming a basis for a broad range of nonlinear optics 
phenomena. Three waves interaction is the simplest representative of such a class of phenomena. 
When frequencies of two waves from an interacting wave triad are equal, the frequency of the 
third wave is two times higher. This is known as second harmonic generation (SHG). Stimulated 
Raman scattering is another example of three- wave interaction. In this case one of the waves is of 
acoustic nature and the incident and scattered waves are electromagnetic. Four wave interaction 
is a more complicated process of parametric interactions. Third harmonic generation (THG) is 
a typical representative of these processes. 

3.1 Three-wave parametric interaction 

Let as assume that the non-linear characteristics of a medium to be described by non-linear 
susceptibility of second order x^^^. Such a medium is commonly known as a quadratic non- 
linear medium. We consider waves with carrier frequencies wi and uj2 propagating along the z 
axis. The polarization of such a medium is a quadratic function of the electrical field, therefore 
waves with frequencies LOg — uji±uj2, 2^1, and 2lu2 must be generated in such a medium. These 
waves, when their amplitudes increase sufficiently to be involved in the process of nonlinear 
interaction, can generate new waves with the frequencies 2uji ± 0^2, i^i i 2w2, 4aji, and 4^2 so 
on. However in a dispersive medium all these processes are not equally efficient. There is an 
important phase matching condition, which selects and emphases a certain type of interaction 
of three waves, leaving all others unaffected. Sometimes such a phase matching occurs for waves 
propagating in the same direction - the case of coUinear parametric interaction. In this case the 
distance at which the interaction of waves is taking place can be made sufficiently long and, 
consequently, the effective frequency transformation will take place. 
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The equations describing three wave interaction {uj^ =0^1+ ^2) in a x^^^-medium in the 
slowly varying envelope and phase approximation, taking into account the group-velocity dis- 
persion, can be written in following form 6J: 

* jj^PNL{iui)expi-ihz) (41) 

C fC]^ 

. |^^L^p^^^2 exp -ifc22 42 

. 27ra;|/i(ijJ3) - 



oz 


1 










h — 

Vl 


dtj 


2 




oz 


1 


a^ 


]S2- 


D2 




h- — 


dtj 


2 




i ^ 


1 


d\ 




D3 




oz 


h — 

V3 


dtj 


1^3- 


2 





PNL{L^3)exp{-ik3z). (43) 



Here /cj is defined as /c| = {ujj /c)'^e{ujj)^{u!j); kj is sign of square root of fc|; and 

pNLit^i) = x*-^''(^i;t^3, -^2)^3^2 exp[iz(fc3 - /C2)], 
Pnl{uj2) = x''^H^2;^3, -^1)^3^1* exp[zz(fc3 - /Ci)], 
PNLiuJa) = x*^''('^3;wi,W2)fif2exp[iz(fci + fc2)]. 

If uji — L02, t-^s — 2wi we have second harmonic generation. If uji — LUg, lu2 — uji, and 0^3 = ujp 
{ujp — ujs +uji), than we are dealing with parametric amplification phenomena. The signal wave 
is amplified by taking energy from the pump wave with the help of the idler wave, playing the 
role of a mediator in this energy transfer. The idler wave here corresponds to £2 and the pump 
wave corresponds to £3. 

The zero phase mismatch condition Ak ~ kp — kg — ki means that the vector of the pump 
wave kp is equal to the vector k^ -I- k^ . The configuration of these vectors defines the energy 
flow directions. If the nonlinear medium is characterized by a positive refraction index, then 
the vector kp and the Poynting vectors associated with the interacting waves have the same 
orientation. Opposite directionality takes place if the nonlinear medium is characterized by 
negative refraction index [9T36|l37|38j . Currently the negative refractive property in the optical 
domain has been realized by use of simultaneous resonance for electric and magnetic field 
components in metallic nanostructures embedded to a host medium. The resonance frequencies 
in these materials are close to each other and they divide the frequency range into two domains. 
The medium responds to external waves as a negative refractive index material if the wave 
carrier frequency is below both resonance frequencies. The medium responds to external waves 
as a positive refractive index material if the wave carrier frequency is above both resonance 
frequencies. In the case of a three wave interaction, wave vectors corresponding to the waves 
with frequencies from the frequency region of negative refraction index are directed according 
to the zero phase mismatch condition, but the associated Poynting vectors are have opposite 
orientation. 

Let us consider parametric amplification. It is convenient to introduce normalized variables 
£1 = Aq-JiCs , £2 = ^07262, £3 = Ao72ei,where jj = 2TruJix^'^'^ {u;i)c~^ ^y ij.{ujj)/e{ujj). The 
system of equations describing the parametric amplification can be rewritten as follows: 



d_ i_d_ 

dz Vl dt 
d_ ^_d_ 
dz V2 dt 
d_ ^ ]_d_ 

dz vs dt 





d^ei 


2 


df^ 


Di 


d'e2 


2 


dt^ 


Di 


d^e. 


2 


dt^ 



ei -t- ^^~Q^ = We2es exp(-iZ\fcz) (44) 

Di 5^62 

62 -I- i— — igeie* exp(-f iZ\/cz) (45) 

igeiBj exp(-f zZ\/cz). (46) 



Here Ak is defined as Ak = kp ~ kg ~ ki and g = (717273)^^^^0 is a coupling constant. The 
value Aq is introduced for appropriate normalization of electric fields of interacting waves. 

In the general case the phase matching condition requires pump and signal waves to prop- 
agate in opposite directions. However, solitary wave solutions of ([15)) are possible if the power 
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of the interacting waves is sufficiently high. Such solitary wave solutions can be considered as a 
bounded state of signal, idler and pump waves propagating in the same direction. This bound 
state can be interpreted as a wave trapping phenomena, where the group velocity of the signal 
wave changes sign due to a nonlinearly induced change in the value of the effective refractive 
index. Note that if the phase mismatch is zero Ak = and the effect of group-velocity disper- 
sion is negligible (short sample) , then equations (|45p reduce to a system of equation integrable 
by use of the inverse scattering transform [32133134] . Therefore, results previously obtained for 
the integrable case of three-wave interaction can also be used to analyze this special case of 
three- wave interaction in LHM. 

Negative refractive index materials based on plasmonic resonance in metallic nanostructures 
have large losses. The feasibility of parametric amplification as a means to compensate for 
dissipative losses in LHM was considered in [39]. It was shown that parametric amplifiers are 
capable of compensating for losses and additionally that they can be used as optical parametric 
oscillators. The opposite directionality of parametrically interacting waves in LHMs is an analog 
of distributed feedback. In conventional materials distributed feedback is realized by adding 
Bragg gratings. 

The case of three-microwave interaction for a strong pump wave and two weak signals 
in metamaterials with nonlinear magnetic response was theoretically studied in [35] , In this 
case nonlinear properties of the material were determined by the embedded LC circuits with 
nonlinear resistance and capacitance. 



3.2 Second Harmonic Generation 

Second harmonic generation is a special case of three wave interaction [loi + loi =0^2) in a 
X*-^-* -medium. The equations of SHG in the case of coUinear mismatch conditions Ak = 2ki — ^2 
follows from (|45)). if Cs is chosen as = ei for fundamental wave (with carrier frequency uj) 
and if ei in ()45p is replaced by 62- In this new notation idler and pump waves correspond to a 
second harmonic: 

- ei,z -I- wf^ei,* -I- i(£)i/2)ei,tt = ige2el exp(-iZ\fcz) (47) 
e2,z + V2^e2,t + i{D2/2)e2,tt = igef e:x.-p{+iAkz) (48) 

Transformations of variables ei — AiqCi, 62 — 2^1062 exp(+i/\fcz) allow us to represent these 
equations in the following standard form 

iei^C + W2)ei,rr-e2et = 0, (49) 
ie2.x + i^^2.T - (/3/2)e2,rr - Ae2 + el/2 = 0. (50) 

Here ^10, C ^tnd t are defined as Aiq = {gL)~^, ( — z/L, t ~ {t + z/vi)/tp. Normaliz- 
ing characteristic parameters here are chosen as follows: L = tp/\Di\ is dispersion length; 
S — Ltp^[vi^ + V2^) is normalized group velocity mismatch; A = AkL is normalized phase 
mismatch; tp is characteristic time, which is not fixed yet; and a — sgn_Di, /? = D2/\Di\. Pa- 
rameter S takes into account the walk-off effect for pump and harmonic pulses that is due to 
the difference of the group velocities' directions for the interaction waves. It should be pointed 
out that in contrast to the case of positive refractive index medium, this parameter can not be 
zero in LHM's. 

The form of the equations (|49|50p is similar to that of the equations describing the evolution 
of quadratic solitons as a bounded complex of fundamental and second harmonic solitary waves 
presented in [40141] . The only difference is in the relative sign of group velocities of the pump 
and second harmonic waves. 

Using the equations (|49l50p and taking into account the boundary condition |ei,2p 
|eio,20p! T ~^ ±00, here Cio and 620 are constants (or |ei^2p — > r ±00), we arrive to the 
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Manley-Rowe relation: 

oo 

J (leal' -(l/2)|ei|2)dr = const. 

— oo 

Note that the sign of the second term in the integrand is negative in contrast to the situation for 
SHG in nonUnear positive index materials. This form of the Manley-Rowe relation reflects the 
fact that the Poynting vectors, i.e., energy fluxes, for the fundamental and the second harmonic 
are antiparallel, while their wave vectors are parallel. 

3.2.1 Continuous wave limit 

The system of equations (|49l50p for continuous waves reduces to following equations [9|37|38j : 

ieix-ele2 = 0, ie2,c - + e^/2 = 0, (51) 

The boundary conditions for a nonlinear plate of a finite width / are as follows: 

At C = |ei|=ao, and at C = ^ |e2| = 0. (52) 

The Manley-Rowe relation in this case reads 

2|e2p - |eip = 2cl = const. 

Solutions of these equations ((5T|) were found in |37|38j . Spatial distributions of the harmonic 
and pump wave are represented by the following expressions: 

|ei(C)| =coy2sec[co(/-C)], |e2(C)| = cq tan [co(/ - C)] • (53) 

The constant cq is defined by the relation aq = coV2 sec [cqCo] that follows from boundary 
conditions (j52p . In the case of quadratic- nonlinear PRI medium the amplitude of the pump wave 
decreases with distance, while the second harmonic wave amplitude increases. The energy of the 
pump wave is transferred into the second harmonic wave and the energy fluxes of both waves 
are aligned. In a quadratic-nonlinear NRI medium the energy fluxes are oriented in opposite 
directions, thus both pump and second harmonic wave amplitudes decrease with distance. 
Second harmonic generation taking into account dissipation has been investigated by numerical 
simulation in the [37j . 

3.2.2 Large-mismatch limit 

During the last few years there has been growing interest in three wave interaction or its special 
case of second harmonic generation when \A\ » 1, (see review j41l). This regime is known as 
the large-mismatch limit, the cascading limit, or the effective Kerr limit. In this limit, equations 
for SHG (|49l50p can be transformed into a nonlinear Schrodinger equation for the fundamental 
wave, and the amplitude of the second harmonic wave is proportional to the squared amplitude 
of fundamental wave. Generally quadratic nonlinearity is an order of magnitude stronger than 
cubic nonlinearity. Therefore the large-mismatch limit case is very useful for study of nonlinear 
cubic phenomena using quadratic nonlinearity. This large-mismatch limit can, in a similar 
way, be considered in the case of a quadratic-nonlinear NRI medium. From (|50p we evaluate 
62 « el/2A, substituting 62 in (|^ results in 

leix + W2)ei,rr - (l/2Z\)|ei|2ei = 0. (54) 

This equation ([51]) has soliton solutions when a — ~l. The type of these soliton solutions 
is controlled by sign of A: these are bright solitons if Z\ > and they are of dark type if 
Z\ < 0. Note that for the frequency of the fundamental (pump) wave the medium acts as left 
handed. On the other hand the governing equation for the pump wave in a large-mismatch 
limit is identical to the conventual NLS equation written for a right handed material. This is 
an indication that the soliton properties of LHM and RHM are similar in the large-mismatch 
limit. 
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3.2.3 Solitary wave solutions 

The solitary wave solutions of (|49l50p . represented by two-frequency pulses, have been consid- 
ered in [42j . To find such solutions it is convenient to introduce real variables for the interacting 
waves 

ei = a exp{itpi) , 62 = 6exp(^(^2)• 
Substitution of these expressions into (|49l50p leads to 

a^(^ + {a/2) {2a_r^i,T + ati^i.Tr) — a5sin^, (55) 
26, (; + 25b - 13 {2b^rf2,T + bip2^Tr) = siu (56) 
O'^i.Q ^ (o'/2) (a.TT — o.^i.tVi,t) = —abcos'P, (57) 
2a {ip2,c + Sip2.t) + P {b^TT + bip2,T^2,T) + 2Ab = a^ cos<P, (58) 

where <P = Lp2 — 2ipi. Assuming that phases are linear functions of the following form tpi — 
K( + Qt, ip2 — 2K( + 2Qt we obtain the mismatch condition = 0. The amplitude equations 
([55)1 and ([5S)) are equivalent if the frequency fi is chosen to he Q — 5/ {a + 2/3). Let us consider 
the solitary wave regime represented by a bound state of two waves with a fixed value of the 
amplitudes' ratio. These waves are defined by functions of a single argument a = a(r — C/^)j 
b = b{T — C/V), with = afi and b = fa. Here / is a constant. Under such an assumption, 
the system of equations (j55p and l|56p is overdetermined (we have two equations for one function 
a — a{y), y = t — C^/V). The compatibility condition for these equation reads: 

f = <j/2p > 0, 

K {3f3n^ - A5Q - 2A) /2{(3 + 2a). 

and the equation for a has the following form 

a^yy - [2aK + f2^)a - 2a fa"^ = 0. 

It follows from the compatibility condition that a and f3 must have the same sign. 

Let us consider the boundary condition a — > 0, da/dy — > at r — > ±00 , which is con- 
sistent with a solitary wave solution propagating on a zero background (quadratic soliton - 
see review [H]). In the case where a = —1 the solitary wave solution of (|55ll58p exists if 
p = {ff' — 2K) > 0. This solution is a bright NRI medium soliton with frequency components 
of the following form 

ajy)^ (3P/4)V2^ ,(^) , (3p/4) ^^^^ 

cosh [^{y - yo)/2] cosh [^{y -yo)/2] 

The group velocity Vg of the bright soliton is fixed by system parameters and defined as 

= {av^'-2Pv^')/{a + 2l3). 

Note, that in the case of a quadratic-nonlinear PRI medium, the group velocity of a bright 
soliton is defined by the formula = (2/3wf^ — crw^^) /(2/? — a). In the case where a = 1 
the solitary wave solution of (|55ll58p exists if p = (i?^ + 2K) > 0. 

If the background is nonzero, then the system of equations (|55ll58p has dark soliton solutions. 
An example of such a dark soliton at cr = —1 and p = {Sl^ — 2K) < is presented below: 

a{y) = (3b|/4)v/2M(2/3±2/Hsech2[/Hy/2]) , (60) 

b{y) = (3b|/4) (2/3 ± 2/H sech2[/Hy/2]) . (61) 

If the ratio of the coupled wave amplitudes is not fixed then the system of equations (I55II58P 
has double-hump solitary waves, which are described by following expression 

a(y) = 3piv/2/3tanh(Vpr2/)sech(^y) , (62) 
b{y) = 3pi sech^ {y/piv) ■ (63) 

where px = (i?^ + 2aK). 
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3.2.4 No group-velocity dispersion limit 

At the current state of the art in nanofabrication technology the interaction distance for para- 
metric three- wave processes is strongly limited by the losses of plasmonic oscillations in metallic 
nanostructures of LHM. The dispersion length is shorter than the characteristic length of losses 
and therefore the dispersion term can be omitted from equations (|49l50p . We also consider the 
case where the length of the LHM sample is shorter than the characteristic length of losses. 
Therefore, the effects of group- velocity dispersion could be omitted as well. The corresponding 
system of equation takes the following form 



The results of SHG computer modeling using these equations were presented in [i^. If Z\ = 0, 
the system of equations (|64l) is integrable using the inverse scattering transform [45] and also can 
be represented in bilinear Hirota form [53] . Note, that in this particular case Hirota's method 
is more convenient to analyze the nonstationary regime of SHG in LHM. 

3.2.5 Non-collinear second harmonic generation 

In this section we consider a process of SHG, where pump waves and the second harmonic wave 
propagate in arbitrary directions. In this case the phase matching condition is determined by 
the beams' orientations. The evolution of slowly varying optical pulse envelopes is described 
by equations, which are similar to (|42p . The spatial derivative in this case must be replaced by 
the corresponding directional derivatives and the argument of the exponential functions in the 
right-hand side must be replaced by ±i(ki + k2 — ks) • r. The phase- matching condition in this 
case reads 



Let us assume that the second harmonic wave propagates along the z-axis - ka vector; both 
wave vectors - ki and k2 - of the pump waves are in the xz- plane; so that ki = ki(rix, —rjz) 
and k2 — ki{~-'qx, —ilz): where rjx = sin9 and yy^ — cos9, 9 is an angle between ki and ks. The 
projection of the vector equation (|65|) onto the z-axis defines the phase-matching angle 9m' 
n(2Lu) = n{uj) cos 0m- Let us assume that the group- velocity dispersion is of no importance. The 
system of equations, describing the SHG, when phase matching condition is satisfied, reads 



ieiX - 626* = 0, ie2x + ^^^2,t - + 6^/2 = 0. 



(64) 



ki + ks - kg = 



(65) 




(66) 



(67) 



(68) 



where the pump wave is presented as 



Ei{t, r) = £[ ' (t, X, z) exp [—iut -I- iki[rixX — rj^z)] + 
+£[ ^ {t, X, z) exp [—iiot — iki{rjxX + rj^z)] , 



and the second harmonic wave is represented as 



E2(t, r) = £2it, X, z) exp [—2iLot -\- ik2z)] . 



System (j67p can be transformed into a system of dimensionless equations 



zei,r + 6362 = 0, «62,i; -I- 6163 = 0, 163,^-1-6^62 = 0, 



(69) 
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where C, and ^ are characteristic co-ordinates defined as 

C = rizViV2ji^/j2AiQ{v2 + rj^viy^ {t - z/rj^vi) , 

r = 2wiU27i\/72^io(«2 + VzViy^ {t ~ z/v2 + {v^;^ + llzV2^)x/'qx) , 

NormaUzed envelopes of the interacting waves are defined by the expressions — ^/Ti^ioei, 

£2 = ^72^1062, S'r^ = VrT^ioea- 

The solutions of these equations could be found either by the inverse scattering method 
[46|47|48] , or by using a combination of simple algebraic manipulations and solution of ordinary 
differential equations 

3.3 Raman scattering process 

Let us consider optical waves propagating in an infinite medium and scattered by optical 
phonons. This process is known as Raman scattering and differs from Mandelschtam-Brillouin 
scattering of optical waves occurs on acoustic phonons. The phonon frequency ujy is much less 
than the carrier frequencies of optical waves. The interaction of optical and vibration modes 
results in a shift of the frequency of the optical mode. This interaction generates two new 
waves (Stokes waves) in addition to the pump wave with carrier frequency uj^. The frequencies 
of these new waves are los — lop — lOv (Stokes wave), and loas ^ + ( anti-Stokes wave). 
In the most cases, the intensity of the anti-Stokes wave is less than that of the Stokes wave. 
Therefore, we consider only incident and Stokes waves. Furthermore, we consider stimulated 
Raman scattering when the frequency up lies in the positive refractive index region and the 
frequency of the Stokes wave lies in negative refractive index region. 

The slowly varying envelopes of the incident and Stokes pulses is governed by the system 
of the reduced Maxwell equations 

/ 9 1 9 \ „ „ 27ra;? uiujp) „ . . , ^ , ^ 

+ j ^ -..8, + ^^^V, eM^Ak.), (70) 

(-^ + -I) - -s£s + ^^-^^^i^VseM-^Ak^l (71) 
\ oz Vg ot J c-^ks 

The linear losses are taken into account by introducing extra terms with the coefficients aL,s- 
We define polarizations Vs.l in the Raman medium in accordance with Placzek's classical 
model where molecules are represented by harmonic oscillators. The polarization of this medium 
is represented as P = nACt{Q)E, where ua is the concentration of the molecules, a{Q) is 
the molecular polarizability, and Q is the vibrational co-ordinate of a molecule determining 
the magnitude of deflection from equilibrium. In the case of small vibrations the molecular 
polarizability is expressed by the first two terms of its Taylor series, i.e., 

a{Q) « ao + {da/dQ)oQ = aa + a^Q- 

Hence, the non-linear polarization is Pnl — nACtoQE. In the slowly varying envelope approx- 
imation Q is represented as follows 

Q = uit, z) exp {—itOvt + ikyz) + u*{t, z) exp {iujvt — ikyz) , 

where u{t, z) is determined by the equation 

u^t — ii^v — ^s)u + <^vU = i {aD/^mujy) S^Es ex.p{iAkz), (72) 

Here m is the effective molecule mass. The coefficient takes into account linear losses. Taking 
into account that the phonon wave number is much smaller than the photon wave numbers, the 
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mismatch value can be taken to be: Ak k, ks — k^. We assume that spectral half-width of the 
ultrashort optical pulse is much smaller than the phonon frequency Aujp <^ Uy and express the 
nonlinear slowly varying amplitudes of polarization as Vl = nAaou* Es and Vs = nAanuSL- 
To develop a theory of Raman scattering we can use the system of equations (I70II72|) . Note 
that in the case we consider the Stokes wave and the laser wave are propagating in opposite 
directions. 

It is convenient to introduce a co- moving system of coordinates = z/L, t = {t — z/Vo)/tp 
with velocity Vq defined by as Vq'^ — (w^^ — Vs^) /2, and normalized functions u = q/(nAaD), 
£l = i/7r^oepexp(iZ\fcz), £3 = ^/isMes, where 

7s = {2tiujsc^^\/ M('^s)/e(ws), IL = {2'!tlopc''^ y/ fi{ujL)/s{ujL), 

The dimensionless system of equations takes the form 

icp.x; + iScp^r — Acp + iFpCp + q*es = 0, 

icsx ~ i^^s^T - iTsCs - qcp = 0, (73) 
iq^r + 'i)q + iTyq + e^e* = 0, 

where 6 = Ltp^{vi^ + V2^)/2, — tp{ujy — ujs), Fp — glL, Fg = agL^Fy — UyL. Therefore 
the system of equations describing Raman scattering is represented as a three wave interaction 
model. 

In the case of strong vibrational damping [Fp ^ 1), it follows from the third equation of 
(f73| that q K, iFy^esCp. Substitution of this formula into the two first equations of (|73|) leads 
to equations describing a Raman amplifier in LHM 

iepx + iSep^r — Acp + iFpCp — iFy^\es\^ep — 0, 
ies^ - iScs^T - iFsCs - iF^^\ep\'^es = 0. 

Note, that choosing Fp = Fg = Q and Fy 1, this system of equations can be used to study 
the phenomena of Raman spike generation in left handed materials. This phenomena is well 
known in conventional right-handed materials |50l51j . 

3.4 Third Harmonic Generation in NRI medium 

Third harmonic generation (THG) is associated with four wave interaction of the type ui + 
wi + tJi — > 0^3 = iiJi . This process is described by the following system of equation 

(fci^; + j ^1 - - ^^l^^-K) exp(-.^i.) (74) 

y^'dl ^V,d-t)^'-^^= ' c-k, "^^^("^^ exp(-^fc3.) (75) 

where 

Pnl{i^i) = x*-^^(t^i;'^3,-wi,-t^i)f3fr^exp[iz(fc3 - 2fci)] -|- 
Pnl{(^2) = x*^^(t^3;'^i,t^i,'^i)fi exp[iz3fci)] -I- 

+X^^\^3;(^3, -W3,CJ3)|f3pf3 + X^^H'^3;'^3, ~UJi, UJi)\£i\'^ £3. 

In contrast to SHG, the process of transformation of one wave to another in this case is accom- 
panied by the two additional effects of self-phase and cross-phase modulation. 
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Let us consider the THG where the frequency of the fundamental wave (pump wave) toi 
is located in the NRI spectral region, and the third harmonic frequency 0^3 is located in the 
positive refractive index spectral region. Let the self-modulation and cross-modulation effects 
be ignored. In this case THG equations after normalization read as 

iei.c + (o-/2)ei,rr - 626*^ = 0, , s 

+ iSe3,r - (/3/2)e3,.r - Acs + e? = 0, ^ "'^ 

where both parameters and normalized variables were introduced along similar lines to this use 
for SHG. 



3.4.1 Third harmonic generation. CW-limit 

First we consider third harmonic generation in the case of continuous waves. The system of 
equations describing THG in the continuous wave limit takes the following form 

ieix - e2ef = 0, ieg,^; - Ae^ -f = 0. (77) 

There are two types of boundary conditions: (a) |ei(C = 0)| = ao,|e3(C — > 00) | = 0, i.e., infinite 
medium, and (b) |ei(C = 0)| = ao,|e3(C = ^)| = , i.e., THG in the plate having width equal to 
I. 

The Manley-Rowe relation following from ((77|) is as follows 

|ei|'- |e3|' = c2 = const. (78) 

Let consider the THG at phase mismatch A = 0. The solution of the system ([77]) for infinite 
medium, can be written as 

|ei(C)l = |e3(C)l = ^==^. (79) 

Hence, in this special case there is complete transformation of the incident wave to a third 
harmonic wave which propagates in the opposite direction. The solution of ()77|) in the more 
realistic case, where the size of the nonlinear dielectric sample I is finite, reads as 

|ei(C)l' = ^, |e3(C)l'-^, e = c2(C-0. (80) 

The equation for the unknown constant Cq follows from the Manley-Rowe relation: 

alii - Pet) = cl. 

Note that in contrast to the case of SHG there is only one solution of this equation which 
has physical meaning. In the case of SHG the corresponding equation reads /aosinco/ = V2cqI 
and in general can have several solution [37l38j . Like in the case of SHG, this Manley-Rowe 
relation follows from the new type of phase-matching condition referred in 37 38] as "backward 
phase-matching" . This type of phase-matching is a fundamental feature of wave interaction in 
NRI medium. 



3.4.2 Solitary wave solutions of third harmonic generation equations 

A solitary wave solution (coupled state of two waves with carrier frequencies loi and 0^3) can 
be found as in the case of SHG. Let us introduce the real variables for interacting waves 
ei = aexp(i(pi) and 63 = 6exp(i(/53), where tpi = KC, + i7r, (p^ = 'iKC, + ?>f2T, a = a{T — C,/V), 
b — b{T — C/V). The system of the equations ([75|) reduces to a single equation for the real 
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amplitude a, if the V ^ ^ a[2, frequency [2 = S/{a + 3/3) , and wave number is determined by 
formula 

Solving this equation we found the field amplitudes 



|ei(T,C)l = . r 1/7/ 7777 TT' |e3(-^,C)l = 



C0sh[pl/2(^ _ ^/V - yo)] ' ' COsh[pl/2(^ _ Q/V - yo)] ' 

where j/q is a constant of integration, — a/ (3 and p = 2<jK + J?-^. For = /9 = — 1 we can 
find Q = -6 /A, K ^ {S^- A)/S and p = (4Z\ - 3(52) ^^q^ rj,-^^ velocity of this sohton is defined 
as T — (/V = {t — z/Vs)/tp and is equal to Vs — 4i;iU3/(z;i — 3u3). 

Note that the system of equations ([76)1 has periodic cnoidal solutions. These solutions de- 
scribe nonlinear periodic waves in a medium with cubic nonlinearity. 

In the large- mismatch limit of THG, the system of equations (j76p transforms into the quintic 
nonlinear Schrodinger equation: 

zei,c + {<j/2)ei.rr ~ {l/2A)\ei\^ei = 0. 
This equation also has solitary pulse like and periodic cnoidal wave solutions. 



3.4.3 Phase interaction of two waves 

There are situations where the process of conversion of waves with different frequencies can be 
ignored. For example, such a situation occurs when the phase difference of the two interacting 
waves is rapidly changing in space: Ak x Ld S> 1, here Ld is the dispersion length. Let us 
consider propagation of two waves in a medium with cubic nonlinearity, where the frequency of 
one wave is in the NRI region, and the frequency of the other wave is in the PRI region. The 
corresponding equations in dimensionless form, which follow from (j74l75[) . read as 

iei.(; (cr/2)ei,rr - (|ei|2 -t-/i|e3|2) ei ==0, , . 

ies.c + iSe^^r - (/3/2)e3,rr + (lesp + /ijeip) 63 = 0, ' 

This system represents a generalization of the Manakov's equations [52 . The solitary wave 
solutions of these equations can be found by standard methods. 



4 Periodic structures 

In recent years, much attention has been drawn to various artificial structures, which are one- 
two- or three-dimensional periodic dielectric (or metallic) structures, frequently referred to as 
photonic crystals. The spectrum of electromagnetic waves in these structures possess bandgaps 
(the Bragg gaps) f53^ . Sometimes the photonic crystals are called photonic band-gap materials 
(PEG). In addition to a conventional PEG, a one-dimensional periodic structure consisting of 
alternating layers of PRI and NRI medium represents a new type of photonic band-gap material, 
a so-called zero refractive index PEG [54j. There we consider only the resonant Eragg grating 
and nonlinear optical waveguide array as simple examples of nonlinear periodic medium. 



4.1 Resonant Bragg grating 

In the simplest case the resonant Bragg grating |55I56I57I58I59| consists of a linear homogeneous 
dielectric medium containing an array of thin films with resonant atoms or molecules. The 
distance between successive films is a, and the thickness of the film If is much less than the 
wavelength of the electromagnetic wave propagating through such a structure. The interaction 
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of ultra-short pulses and films embedded with two-level atoms have been studied by Mantsyzov 
et al. [55156] in the framework of the two- wave reduced Maxwell-Bloch model and by Kozhekin 
[57l58j . This work demonstrated the existence of the 27r-pulse of self-induced transparency in 
such structures. It was also found [60j that bright as well as dark solitons can exist in the 
prohibited spectral gap, and that bright solitons can have arbitrary pulse area. If the density 
of two-level atoms is very high, then the near-dipole-dipole interaction is noticeable and should 
be accounted for in the mathematical model. The effect of dipole-dipole interaction on the 
existence of gap solitons in a resonant Bragg grating was studied in Ref. [6l]. The unusual 
solution known as a zoomeron was discovered and investigated recently in the context of the 
resonant Bragg grating [62 . A zoomeron is a localized pulse similar to an optical soliton, except 
that its velocity oscillates about some mean value. 

Recent advances in nanofabrication have allowed the creation of nanocomposite materials 
which have the ability to sustain nonlinear plasmonic oscillations. These materials have metallic 
nanoparticles embedded in them ^63 64 65 . In [66 a dielectric into which thin films containing 
metallic nanoparticles have been inserted was considered. These thin films are spaced periodi- 
cally along the length of the dielectric so that the Bragg prohibited spectral gap is centered at 
the plasmonic resonance frequency of the nanoparticles. A set of discrete equations describes 
the propagation of light through a medium of alternating layers of linear dielectric and thin 
film. However, the physical parameters of the system make the slowly- varying envelope approx- 
imation appropriate. Following reviews |67|68|69] . we employ such an approximation to derive 
envelope equations for two counter-propagating electric fields which are coupled through the 
medium polarization. The electric field can be described by 

E{x,t) = {Aix,t)e"^°^ + Bix,t)e-"^"^} e-''^"* . 

Here ujq is the carrier frequency, and qo = ujQy/e/c is the wavevector in an optical medium with 
electric permittivity e. The long- wave and slowly varying envelope of the pulses approximations 
result in the system of coupled wave equations. 

« (Ax + V'AO + ^QoA = -(27rwo/cV^)^, (82) 
* - - ^<lo>3 = +(27rwo/cVi)^, (83) 

where Aqo = qo — 27r/a is the mismatch between the carrier wavenumber and the Bragg reso- 
nant wavenumber. Description of the slowly varying polarization is based on the anharmonic 
oscillator model ([55]) : 

V,t - i{L^o - LOd)V-{3igp/2LUo)\P\''V = iiLul/8TrLOo)£f, 

where Smt is the electric field interacting with metallic nanoparticles. In our case we have 
£int = A + B. Note, that if anharmonicity of oscillators is neglected {gp = 0) then the system 
of equations transforms into the well known Lorentz model. In the general case the system of 
resulting equations is the two-wave Maxwell- Duffing equations. Using rescaling 

/, = -{A + B)A^'e-''\ fa^iA- B)A^'er''\ 
q = {ATTUJo/[y^LOpAo])e-"^^V, C = {ujp/2c)x, t = t/to- 

these equations can be represented in dimensionless form: 

/a,CC ^ /a.TT = 2zg_^ 

/s,CC - h,rr = 2iq^r (84) 

iq.r + {A-5)q + lL\q^q= fs. 

Here Iq = l^ju^p is an inverse plasma frequency, is characteristic amplitude of counter- 
propagating fields, \x = (3(7p-\/£/a;oWp)(-\/e[jp/47r[Jo)^^o i^ ^ dimensionless coefficient of anhar- 
monicity, (5 = 2Z\go(c/wp) is the dimensionless mismatch coefficient, A = 2y/e{ujd — uJo)/u}p 



Will be inserted by the editor 



19 



is dimensionless detuning of a nanoparticle's resonance frequency from the field's carrier fre- 
quency. 

The exact solitary wave solutions of the system (|84)) have been found in ^66] It was shown 
that, in contrast to conventional 27r-pulses, they have nonlinear phase. The stability of these 
solutions is sensitive to the phase perturbations. It was also demonstrated that the outcomes 
of pulse collisions are highly dependent on relative phase. 

4.2 Nonlinear optical waveguide array 

A simple example of a waveguide structure is the directional coupler. The directional coupler is 
represented by two parallel waveguides. The separation distance between these waveguides so 
small that guided light is leaking from one waveguide to another. Sometimes these waveguides 
are described as tunnel-connected [8 . If the sign of the refractive index in both waveguides 
is positive, then a beam launched into one of the waveguides produces waves in both waveg- 
uides propagating in the same directions. If the index of refraction of one of the waveguides 
is negative, then a beam launched into one of the waveguides produces two beams in these 
waveguides propagating in opposite directions '70]. The theory of this anti- directional coupler 
is similar to the description of SHG under the undepleted pump approximation. Hence, the 
anti-directional coupler acts as distributed mirror [5]. If the waveguides are made from non- 
linear dielectric material or if they are embedded into a non-linear medium, then the coupling 
performance, in addition to the coupler's geometry and material property, depend on the wave 
input power. If input power is above a certain threshold the antidirectional coupling property 
can be transformed to become unidirectional. 

Nonlinear optical waveguide arrays (NOWA) are a natural generalization of nonlinear cou- 
plers. NOWA with a positive refractive index have many useful applications and are well studied 
in the literature (see for example |71|73|72] ). Materials with negative index of refraction offer 
new types of array structures which can be combined in three groups. The first group corresponds 
to an array of alternating waveguides with different signs of refractive index. The second group 
corresponds to the situation when one planar array containing identical waveguides with pos- 
itive sign of refractive index is attached to another planar array of waveguides with negative 
sign of refractive index. The third group corresponds to a situation where an array of identical 
waveguides contains a single waveguide with opposite sign. 

To develop a theory of such waveguides it is natural to utilize existing formalism based on 
nearest-neighbor coupling. The electric field of an optical wave propagating in NOWA in the 
positive z direction can be represented as follows 



The mode function for a particular m-th mode of channel J is denoted by {x,y), where 
J = 0, ±1, ±2, . . ., and A'd^ (2, t) is a slowly varying envelope of the electric field corresponding 
to this mode. Parameters (3m^ are propagation constants. Omitting the details we can write 
the general equations which are governed by normahzed envelopes = Am\z,t)/Ao : 



Here we assume that the group velocities and the propagation constants are equal for all 
channels and only resonant coupled modes are essential for our further consideration (crj = (jgm 
and Vg = Vgm)- The coefficient K12 is the coupling constant between neighboring waveguides, 
and Xef f is the effective non- linear susceptibility. The amplitude ^0 is chosen to normalize the 
nonlinear coefficient of self-interaction to one. 



J— + OC 



Eix, y,z,t)= J2 ^)'^™ ^ y) {-i^ot + if3l;Pz^ 



J= — oo m 




= 0. 



(85) 
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Using equations ([85]) and taking into account that (fcj ~ +1) and (kj — —1) corresponds 
to PRI and NRI cases respectively, we can present mathematical models describing the three 
groups of arrays listed above. 

4.2.1 Alternated nonlinear optical waveguide array 

Let us suppose that waveguides marked a,s J — 2n and J = 2n + 1 have PRI and NRI proper- 
ties respectively. The system of equations describing electromagnetic wave propagation in this 
structure reads 

iq2n,z + iVg^q2n,t - 'J2nq2n,tt + K12 {q2n-l + 92«+l) + \q2n? q2n = 0, 

«92n+l,2 — iv~^q2n+l,t + (^2n+iq2n+l,tt — K12 {q2n + q2n+2) ^ |<Z2n+l P<Z2n+l = 0. 

In the linear limit these equations correspond to the model of a harmonic lattice with alternated 
sign for the coupling. In this case the spectrum of the linear waves consist of two dispersion 
branches. This means the alternated linear optical waveguide array acts as a gap medium. 
Therefore, in the nonlinear case we can expect the existence the solitary waves. 

4.2.2 Interface of the nonlinear optical waveguide arrays 

Let suppose that waveguides marked asn > and n < have PRI and NRI properties respec- 
tively. In this case propagation of electromagnetic waves is governed by the following system of 
equations 

iqn,z + iVg^qn,t ~ CFnqnM + Ki2 (g„-l + Qn+l) + |gnp(7n =0, n > 0, 
iqn.z - iVg^qn,t + <7nqn,tt - Ki2 {qn-1 + qn+l) - knPgri =0, n < 0. 

In the linear approximation, due to the opposite directionality of phase and group velocities 
in NRI materials, one can expect existence of surface waves with vortex structures. Each type 
of nonlinear infinite array supports solitary wave solutions. One can expect that the boundary 
effect at the interface of two types of arrays will be revealed through pinning of solitary wave 
solutions in a way analogous to the behavior of light near defects in a photonic crystal. 

4.2.3 Nonlinear optical waveguide array with a defect 

Let us assume that the PRI type array of waveguides contains one waveguide of NRI type. This 
waveguide, marked n — 0, can be viewed as a defect in a one dimensional lattice. Propagation 
of an electromagnetic wave in such an array is governed by the equation: 

iqn.z + iVg^qn.t - CTnqnM + K12 (g„-l + qn+l) + \qn\'^qn =0, n > 0, 

iqo,z - iVg^qo.t + croqoM - K12 {q-1 + q+i) - \qo\^qo = 0, n = 0, 

iqn.z + iVg^qn.t - CTnqnM + K12 {qn-1 + qn+l) + l^nPgn =0, n < 0. 

Presence of a "defect" could be a reason for a pinning phenomena, where the electromagnetic 
wave is trapped in this waveguide. One can expect that increase of the localized field energy 
will destroy localization. 

5 Conclusion 

We considered new features of classical phenomena of nonlinear optics in optical metamaterials 
including materials with negative refractive index. In addition to already obtained results in 
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this field, we discussed a number of problems of potential importance for the nonlinear optics 

of mctamaterials. 

Optical transparency of metamaterials is a fundamental assumption in our considerations. 
The current state of the art in nanofabrication technology is capable of delivering only thin 
film metamaterials with a relatively large level of losses. We expect that further development 
in nanofabrication will overcome problems of metamaterial losses. Our expectations are based 
on the history of optical fiber technology development in the second part of the last century. 



Acknowledgment 

We would like to express our gratitude to J-G. Caputo, R. Indik, N.M. Litchinitser, B.I. Mantsy- 
zov, and A. A. Zabolotskii for enlightening helpful discussions. A.I. Maimistov thanks the De- 
partment of Mathematics of University of Arizona and Laboratoire de Matliematiques, INSA de 
Rouen for the support and hospitality while he looked into the problems under consideration. 
This study was supported in part by the Russian Foundation for Basic Research (grant no. 
06-02-16406), by NSF (grant DMS-0509589), ARO-MURI award 50342-PH-MUR and State of 
Arizona (Proposition 301). 



References 

1. R. Shelby, D. R. Smith and S. Schultz, Science, 292, (2001) 77. 

2. L.I. Mandol'shtam, Zh. Exp. Thcor. Fiz., 15, (1945) 475 

3. V.G. Vesclago, Sov. Pliys. Uspckhi 10, (1968) 509-514 

4. J.B. Pendry, Phys. Rev. Lett., 85,(2000) 3966-3969 

5. J.B. Pendry, Phys. Rev. Lett. ,85, (2000) 3966-3969 

6. Shen, Y.R. The principles of non-linear optics (John Wiley Sons, New York, Chicester, Bris- 
bane, Toronto, Singapore, 1984). 

7. R. Ruppin, Phys.Lett. A 277, (2000) 61-64 

8. Guided-Wave Optoelectromcs, ed. T. Tamir (Springer- Verlag, Heidelberg, 1988). 

9. V. M. Agranovich, Y.R. Shen, R.H. Baughman, A. A. Zakhidov, Phys. Rev. B 69, (2004) 165112 

10. S.A. Darmanyan, M. Neviere, and A. A. Zakhidov, Opt.Commun. 225, (2003) 233-240 

11. S.A. Darmanyan, M. Novioro, and A.A. Zakhidov, Phys.Rev. E 72, (2005) 036615 

12. R. Ruppin, J. Phys.: Condens. Matter 13, (2001) 1811-1818. 

13. I. V. Shadrivov, A. A. Sukhorukov, and Yu. S. Kivshar, Phys.Rev. E 67, (2003) 057602. 

14. I.V. Shadrivov, Photonics Nanostruct. Fundam. Appl. 2, (2004) 175-180. 

15. G. D'Aguanno, N. Mattiucci, M. Scalora, and M. J. Bloemer, Phys.Rev. E 71, (2005) 046603. 

16. J. Schollong, C. Monzon, P. F. Loschialpo, D.W. Forester, and L. N. Medgyesi-Mitschang, Phys.Rev. 
E 70, (2004) 066606. 

17. De-Kui Qing, Gang Chen, Opt.Lett. 29, (2004) 872-874. 

18. I. V. Shadrivov, A. A. Zharov, and Yu. S. Kivshar, Appl.Phys.Lett. 83, (2003) 2713-2715. 

19. A.A. Zharov, N.A. Zharova, I.V. Shadrivov, Yu.S. Kivshar, Appl.Phys.Lett. 87, (2005) 091104. 

20. V.I. Rupasov, and V.I. Yudson, Sov.J.Quantum Electron. 12, (1982) 415-419. 

21. V.I. Rupasov, and V.l. Yudson, Sov. Phys. JETP 66, (1987), 282-285. 

22. A.M. Basharov, Sov.Phys. JETP 67, (1988) 1741-1744. 

23. E.Vanagas, and A.I. Maimistov, Opt. Spektrosk. 84, (1998) 301-306. 

24. A.M. Basharov, A.I. Maimistov, and E.A. Manykin, Sov.Phys. JETP 70, (1990) 864-871. 

25. I. R.Gabitov, R.A. Indik, N.M. Litchinitser, A.I. Maimistov, V.M. Shalaev, J.E. Soneson, J. Opt. 
Soc. Am. B 23, (2006) 535-542 

26. H. M. Gibbs, Optical Bistability, (Academic Press, Orlando) 1985. 

27. Y. Ben-Aryeh, CM. Bowdcn, J.C. Englund, Phys.Rev. A 34 , (1986) 3917-3926. 

28. Th. Peschel, P. Dannbcrg, U. Langbein, F. Lederer, Opt. Soc. Am. B 5, (1988) 29-36. 

29. C.Z. Koopko, Phys. Stat. Solidi (b) 141, (1987) K139-K142. 

30. N.M. Litchinitser, I.R. Gabitov, A.I. Maimistov, and V.M. Shalaev, Opt.Lett. 32, (2007) 151-153. 

31. I.V. Shadrivov, Photonics Nanostruct. Fundam. Appl. 2, (2004) 175-180. 



22 



Will be inserted by the editor 



32. V.E. Zakharov, and S.V. Manakov, Zh.Eksp.Teor.Fiz. 69, (1975) 1654-1673. 

33. D.J. Kaup, Stud. Appl. Math. 55, (1976) 9-44. 

34. M.J. Ablowitz, and H. Segur, Solitons and the Inverse Scattering Transform (SIAM, Philadelphia, 
1981). 

35. M. Lapine and M. Gorkunov, Phys.Rev. E 70, (2004) 066601. 

36. I.V. Shadrivov , A. A. Zharov, Yu.S. Kivshar, Second-harmonic generation in nonlinear left-handed 
metamaterials. arXiv: physics/0506092 J.Opt.Soc.Amer. B. 23, (2006) 529-534. 

37. A.K. Popov, V.V. Slabko, and V.M. Shalaev, Laser Phys. Lett. 3, (2006) 293-297. 

38. A. K. Popov and V. M. Shalaev, Appl. Phys. B 84, (2006) 131. 

39. A.K. Popov, V.M. Shalaev, Opt. Lett. 31, (2006) 2169-2171. 

40. B. Malomed, Solitons in optical media with quadratic nonlinearity in Nonlinear Science at the 
Dawn of the 21st Cen-tury, ed. by P.L. Christiansen, M.P. Sorensen, and A.C. Scott, (Springer, 
2000) pp. 247-262. 

41. A.V. Buryak, P. Di Trapani, D.V. Skryabin and S. Trillo, Optical solitons due to quadratic nonlin- 
earities: from basic physics to futuristic applications, Physics Reports 370, (2002) 63-235. 

42. A. I Maimistov, L R. Gabitov, E.V. Kazantseva, Optics and Spectroscopy 102, (2007) 90-97. 

43. M. Scalora, G. D'Aguanno, M. Bloemer, M. Centini, D. de Ceglia, N. Mattiucci, and Yu. S. Kivshar, 
Optics Express 14, (2006) 4746-4756. 

44. R. Hirota, J.Phys.Soc.Japan 46, (1979) 1927-1928. 

45. D.J.Kaup, Stud.Appl.Math. 59, (1978) 25-35. 

46. D.J. Kaup, Physica D 1, (1980) 45-67. 

47. D.J. Kaup, Stud.Appl.Math. 62, (1980) 75-83. 

48. D.J. Kaup, Physica D 3, (1981) 374-395. 

49. A.D.D. Craik, Proc.Roy.Soc. (London) A363, (1978) 256-269. 

50. K. Driihl, R.G. Wenzel, J.L. Carlsten, PRL 51, (1983) 1171 

51. R.G. Wenzel, J.L. Carlsten, K. Driihl, J. Stat. Phys., 39 (1985) 621 

52. S.V. Manakov Sov.Phys. JETP 38 (1974 ), 248-256 

53. J. D. Joannopoulos, R.D. Meade, J.N. Winn, Photonic crystals (Princeton University Press, Prince- 
ton 1995). 

54. J. Li, L. Zhou, C.T. Chan, and P. Sheng, Phys.Rev.Lett. 90 (2003), 083901. 

55. B.L Mantsyzov, D.O. Gamzaev, Optics and spectroscopy 63 (1987), 200-202 

56. B.L Mantsyzov, Phys.Rev. A 51, (1995) 4939-4943. 

57. A. Kozhekin, G. Kurizki, Phys.Rev.Lett. 74(1995) 5020-5023. 

58. A. Kozhekin, G. Kurizki, B.Malomed, Phys.Rev. Lett. 81 (1998), 3647-3650. 

59. G. Kurizki, A.E. Kozhekin, T. Opatrny, B. A. Malomed, Progr. Optics. 42, (E. Wolf, editor: North 
Holland, Amsterdam, 2001) pp. 93-146 

60. T. Opatrny, B.A. Malomed, G. Kurizki, Phys. Rev.E 60 (1999), 6137-6149. 

61. J. Cheng, J. Zhou, Phys.Rev. E 66 (2002), 036606. 

62. B.L Mantsyzov, JETP letters 82 (2005), 253-258. 

63. S. G. Rautian, JETP 85 (1997) 451-461. 

64. V. P. Drachev, A. K. Buin, H. Nakotte, and V. M. Shalaev, Nano Lett. 4 (2004), 1535-1539. 

65. F. Hache, D. Ricard, and C. Flytzanis, J. Opt. Soc. Am. B 3 (1986), 1647-1655. 

66. A.L Maimistov, LR. Gabitov, A.O. Korotkevich, Quamtum Electronics (2007) in press. 

67. C. M. De Sterke, J. E. Sipe, Prog. Opt. 33 (1994), 203 . 

68. R. Kashyap, Fiber Bragg Gratings (San Diego: Academic Press, 1999). 

69. G. Agrawal, Y. S. Kivshar, Optical Solitons: From Fibers to Photonic Crystals (Amsterdam: Aca- 
demic Press. 2003). 

70. A. Alu and N. Engheta, in Negative-refraction Metamaterials: Fundamental Principles and Appli- 
cations, eds. G. V. Eleftheriades and K. G.Balmain (Wiley, New York, 2005). 

71. D.N. Christodoulides, and R.I. Joseph, Opt. Letts. 13 (1988), 794-796. 

72. S. Darmanyan, I. Relke, and F. Lederer, Phys.Rev. E 55 (1997), 7662-7668 

73. C. Schmidt-Hattenberger, U. Trutschel, and F. Lederer, Opt.Letts. 16 (1991), 294-296. 



